Scaling at chiral quantum critical points in two dimensions 
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We study the locahzation properties of electrons moving on two-dimensional bi-partite lattices in 
the presence of disorder. The models investigated exhibit a chiral symmetry and belong to the chiral 
orthogonal (chO), chiral symplectic (chS) or chiral unitary (chU) symmetry class. The disorder is 
introduced via real random hopping terms for chO and chS, while complex random phases generate 
the disorder in the chiral unitary model. In the latter case an additional spatially constant, per- 
pendicular magnetic field is also applied. Using a transfer-matrix-method, we numerically calculate 
the smallest Lyapunov exponents that are related to the localization length of the electronic eigen- 
states. From a finite-size scaling analysis, the logarithmic divergence of the localization length at 
the quantum critical point at J5 = is obtained. We always find for the critical exponent k, which 
governs the energy dependence of the divergence, a value close to 2/3. This result differs from the 
exponent k = 1/2 found previously for a chiral unitary model in the absence of a constant magnetic 
field. We argue that a strong constant magnetic field changes the exponent n within the chiral 
unitary symmetry class by effectively restoring particle-hole symmetry even though a magnetic field 
induced transition from the ballistic to the diffusive regime cannot be fully excluded. 

PACS numbers: 71.23.An,73.20.Fz,72.80.Vp 



I. INTRODUCTION 

The properties of electrons moving in disordered two- 
dimensional (2D) lattices featuring special symmetries 
represent a long-standing challenge for theories of An- 
derson localizationji It is known that the basic mech- 
anism for electron localization, i.e., quantum interfer- 
ence of an electron moving along different paths in disor- 
dered media, depends on the symmetry and topology of 
the particular modelj^i^ Recently, the discovery of topo- 
logical insulators^ and the utilization of graphene^*^ — a 
two-dimensional honeycomb lattice of carbon atoms — 
has been producing interest in this topic anew. In the 
presence of disorder that preserves the sublattice sym- 
metry of the underlying two equivalent interpenetrating 
triangular sublattices, the energy eigenvalues still come 
in pairs ie^ around the Dirac point at = 0. The phe- 
nomenon of a quantum critical point at zero energy is 
not linked to graphcne's hexagonal lattice but can arise 
also in other bi-partite lattices like the simple square lat- 
tice. Disordered systems possessing this property belong 
to the chiral symmetry classesi^ One distinguishes three 
chiral Gaussian random matrix ensembles: the chiral or- 
thogonal (ChO), the chiral unitary (ChU), and the chi- 
ral symplectic (ChS) one, which differ from each other 
by the presence (ChO) or absence (ChU) of time reversal 
symmetry. The chS ensemble has time reversal symmetry 
but spin-rotational symmetry is lacking, due to spin-orbit 
interactions in most cases. An additional particle-hole 
symmetry is obeyed in the chO and chS classes. 

In disordered 2D chiral systems, all electronic states 
are localized with an energy dependent localization 
length that diverges at the band center E = 0.—"— The 
predicted divergence of the disorder averaged density 
of states (DOS) obtained from non-linear sigma-model 



field theoriesi^^— was apparently observed in numerical 
studieSfi^ although the precise value of the critical ex- 
ponent remained unclear. On the contrary, in numer- 
ical investigations^^ of Motrunich et al. for chiral or- 
thogonal systems (time-reversal symmetry is conserved) , 
the exponent k of the asymptotic energy dependence of 
the disorder averaged DOS, p{E) ^ exp(— c| \nE\'^) 
was found to be compatible with k = 2/3, where c is 
a nonuniversal constant. This result is at variance with 
the predicted universal exponent k = 1/2^^^"— that 
is agreed to be correct to all orders of perturbation the- 
ory. At first, this difference was believed to be due to 
the strong disorder limit applied in Ref. (|20l ). The same 
result K = 2/3 was obtained, however, also in the limit of 
weak disorder— for a model of non-interacting particles 
moving on a square lattice via random nearest neighbor 
hopping in the presence of a 7r-fiux phasei^ Very recently, 
it has been pointed out that non-perturbative effects re- 
lated to topologically non-trivial excitations have to be 
taken into account in field-theoretical investigations of 
quantum phase transitions in disordered systems belong- 
ing to the chiral symmetry classesi^ 

Besides the divergence of the density of states dis- 
cussed above, a corresponding divergence of the localiza- 
tion length with k = 1/2, also predicted by theory^! to 
occur in a very narrow energy interval close to E — 0, has 
been successfully observed for a chU model in a numerical 
study of the smallest Lyapunov exponent only recently?^ 
despite several previous attemptsi^"— Scaling properties 
of 2D chiral systems can be studied numerically by the 
analysis of the length and disorder or energy dependence 
of Lyapunov exponents defined for quasi-one-dimensional 
(qlD) systems. In the numerical analysis, strips of fi- 
nite width L and length ^ L are calculated by the 
transfer matrix method. The quantity of interest is the 
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smallest eigenvalue xi of the matrix ln(T'''T), where T 
is the transfer matrix. Following the finite size scaling 
method, from which the localization length of the infi- 
nite 2D system can be inferred by studying qlD samples, 
we assume that the smallest positive Lyapunov exponent 
7 = limL_^oo xiL/Lz is a function of the ratio InL/ In^, 
where ^{E) is the correlation length. 

Using a bricklayer lattice model with uncorrelated ran- 
dom-magnetic-flux (RMF) disorder;^ distributed uni- 
formly with zero mean, a scaling of the dimensionless 
quantity zi = xiL/L^ was observed;^ 



ziiE,L) = Fi 



InL 



(1) 



with the scaling function Fi depending only on the ratio 
(In L) / (In ) . This model belongs to the chiral unitary 
symmetry class (broken time- reversal symmetry). The 
same exponent was found also for a similar model defined 
on the square latticei^ This result confirms theoretical 
predictions that for \E\ <C Eq the localization length is 
^{E) = ^0 exp[^(ln(£'o/|i?|))''] with a universal exponent 
At = l/2ii^i2i However, the recent suggestion^! that non- 
perturbative effects have to be taken into account applies 
also to these investigations. 

It is the aim of our work to improve the understand- 
ing of critical properties at chiral quantum critical points, 
and to find possible reasons for the differing critical expo- 
nents obtained in the energy dependence of the DOS and 
in the divergence of the localization length. Therefore, 
we present in this paper further comprehensive numer- 
ical results obtained for different disorder strengths in 
model systems exhibiting chiral orthogonal, chiral sym- 
plcctic, or chiral unitary symmetry in the presence of a 
very strong perpendicular magnetic field. The latter sit- 
uation is quite intriguing because it is well known that in 
a 'normal' two-dimensional (diagonal) disordered system 
with broken time-reversal symmetry (unitary) all elec- 
tronic states are localized. However, this situation con- 
verts into the quantum Hall case, where current carry- 
ing critical states appear in the disorder broadened Lan- 
dau bands, if a strong normal magnetic field is applied. 
Therefore, the question arises whether in the case of chi- 
ral unitary symmetry the properties in the vicinity of the 
chiral critical point at E = are also influenced by the 
presence of a strong magnetic field. 

We find in all three cases mentioned above a diver- 
gence in the energy dependence of the localization length 
at ii^ = with a critical exponent k ~ 2/3 similar to the 
exponent found in Refs. (j20ll22[ ) for the diverging disor- 
der averaged DOS. This value is, however, in contrast 
to our previous result for the chiral unitary case without 
a constant magnetic field, where k ~ 1/2 was obtained. 
We argue that a broken particle-hole symmetry may be 
responsible for this difference. 



II. MODELS 

We consider two dimensional (2D) tight binding mod- 
els defined on bi-partite square or bricklayer lattices with 
periodic boundary conditions applied in both directions. 
The latter model has recently been shown to be well 
suited for describing the electronic properties of non- 
interacting particles on graphene's honeycomb lattice in 
the presence of disorder Also, in contrast to the square 
lattice, the van Hove singularities do not coincide with 
the quantum critical point at E ~ which facilitates 
the observation of the energy dependence of the criti- 
cal divergences. In the absence of diagonal disorder the 
Hamiltonian reads 



n 



(2) 



where cj^ and Cn are the fermionic creation and annihi- 
lation operators at lattice site n, and the sum runs over 
nearest neighbors only. The disorder is incorporated in 
the real random hopping (RRH) terms 
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to exp —e 
to 



(3) 



pointing in the transversal direction and connecting every 
other pair of atoms (those bonds between the zigzag lines 
on a hexagonal lattice) in the bricklayer situation. Here, 
{s} is a set of uncorrelated random numbers with box 
probability distribution, |e| < 1/2. In our investigation, 
the strength of the disorder is varied between W/to = 2 
and W/to = 8. The spectral bandwidth A is determined 
by the sum over the mean of the hopping terms and the 
coordination number Z, A = {Z/2)({tx) + {ty)), with 
tx/to = 1 and 



(ty) = io / de exp 

J-l/2 



W 
— e 

to 



W 2to 



(4) 



A increases by a factor of 6 when the disorder strength 
W is increased from W/to = 4 to W/to = 8. 

The same off-diagonal disorder type with W/to = 4 
was used in the chiral symplectic model defined on an 
ordinary square lattice. We consider a chiral version 
of the two-dimensional Ando model.— The hopping 
terms i„„' are now 2x2 matrices 



where + = 1, s = 1/2, and the disorder in the 
hopping ty is given by Eq. ([3]). 

The chiral unitary model is studied again on a brick- 
layer lattice^i The hopping terms in the transversal di- 
rection are defined as 



ty = toe'^^-y^^-y^ 



(6) 
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where the phases 9x,y;x,y+a = Ox+2a,y;x+2a,y+a " ^^x.y 

are determined by the total flux threading the plaquette 

at (x, y) 



p h 
q e 



(7) 



where p/q {p and q are mutually prime integers) is the 
rational number of magnetic flux quanta h/ e per plaque- 
tte 2a? , and B = ph/{2qea?') is the magnetic flux density 
perpendicular to the two-dimensional lattice. This dif- 
fers from the random flux model of our previous work2^, 
where the constant magnetic field part was absent. The 
random part originates from the magnetic fluxes 4>x,y^ 
which are uniformly distributed —jjl < (f)x,y < //2 with 
zero mean and disorder strength < f /{h/e) < 1. 



III. METHOD 

Wc calculate the two smallest dimensionlcss scaling 
variables zi , Z2 for quasi-one-dimensional samples having 
an even width L in the range 8 < L/a < 160 lattice spac- 
ings a using the transfer-matrix-methodi^ For E = 0, 
z\ =2L/\\s related to the electronic localization length 
\{E = 0, L). The values of zi, Z2 were obtained with the 
relative accuracy ^Juxj z\ ~ a x 10~^, where a = 3 for 
small L, and a increases to 10 for Lja = 128 and larger. 
For zero energy, we checked that the difference 



zi(£:=:0,i)-Z2(£: = 0,i) 



(8) 



is smaller than the numerical uncertainty of our data for 
all L. For the energy dependence, a narrow interval of 
\E\ < 10~®io around the critical point at E = was 
considered. 

While investigating various disorder strengths increas- 
ing from W/ta = 2 to W/to = 8 for the chiral orthog- 
onal model, it turned out that the scaling is more dif- 
ficult to observe for weak disorder. The reason is that 
in this case the smallest scaling variable zi is less than 
1, which means that the related localization length ex- 
ceeds the sample width. Therefore, putting, for instance, 
W/to = 2, we could not observe the insulating regime for 
non-zero energies even on length scales L/a < 160, and 
our data suffer from strong flnite-size effects. Therefore, 
in what follows we discuss only data for stronger disorder 
W/to > 5 with zi > 1. 

Wc flt our numerical data to the polynomial functions 



Zi,2{E, L) = Zi,2{E = 0, L) ± ^ hn 



X 



(9) 



where Zi{E = 0, L) and Z2{E = 0, L) are the numerically 
obtained values related to the two smallest Lyapunov ex- 
ponents at zero energy. The dimensionlcss scaling pa- 
rameter X is defined as 



X 



[\n{L)/AY^ 
[HEo/\E\)Y^- 



(10) 



We assume the following energy dependence of the cor- 
relation length 



Hm/io) = A\\n{Eo/\E\r, 



(11) 



with unknown parameters i?o, ^Oi a-nd critical exponent 

K = X2/Xz. 

The total number of fitting parameters is Np — A + 
n f , where n j determines the highest order of polynomial 
in Eq. ([9|). The additional four parameters are the two 
exponents X2 and X3 , the energy Eo , and the length scale 
^0- The latter was fixed to ^0 = 1 in all our analysis.— In 
the fitting procedure, we minimize the following function 



F = 



1 ^ 
N ^ 



[z^^\E,L)-Z^iE,L)Y^ 



An) 



iE,L) 



cTiiE,L) 
Z2{E,L)Y 



(12) 



In Eq. (H 



a2{E,L) 



and Z2"' are the numerically obtained 



data for the first and the second smallest scaling variable, 
and CTi and 02 are their numerical uncertainties. To es- 
timate the accuracy of the fit, we create a statistical en- 
semble of up to A^stat = 50 sets with initial conditions 
taken as Zi + y^oT x e with random number |e| < 1. 

After fitting the obtained numerical data for zi and Z2 
to the function 0, we change the number of fitting pa- 
rameters (i.e., the order n/ of the polynomial) in order 
to check the reliability of the fit. Although we do not 
expect strong finite size effects, we also vary the smallest 
size of the systems taken into account in the data en- 
sembles. Another check of the stability of the resulting 
exponent is to reduce the energy interval, for instance to 
values E < 10-^%. 



IV. RESULTS 

The results for all models analyzed are summarized 
in Table HI For the chiral orthogonal model, we present 
data for four different values of the disorder strength in 
the range 5 < W/to < 8. These disorders are strong 
enough to assure that the smallest scaling variable is zi > 
1, so that we are in the localized regime for non-zero 
energies. The exponents k are close to the value 2/3. 
Small deviations are seen for weak disorder W = 5. 

For the chiral symplectic system, wc analyzed only one 
value of the disorder strength, W/to ~ 4. The obtained 
value for the critical exponent is compatible with k = 
2/3 and agrees well with the one found for the chiral 
orthogonal symmetry. 

For chiral unitary systems, we consider a model with a 
constant magnetic flux density B = (1/12) h/(2ea'^) and 
RMF disorder strength f/{h/e) = 0.5 (labeled chU[l/12, 
0.5]), and also a model where the constant magnetic field 
B = (1/16) /i/(2ea^) is combined with real hopping terms 
ty as given by Eq. (jS]) (model chU[l/16, 5.0]). Both sys- 
tems exhibit a critical exponent consistent with a value 
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TABLE I: The critical parameters for models with RRH 
or RMF disorder with a strong constant magnetic field as 
obtained by the fitting procedure described in Section IIIII 
In the ensembles marked by *, a restricted energy interval 
\E\ < 10"^^io was used instead of \E\ < 10"*to- The labeling 
of the chiral orthogonal chO[ ], symplectic chS[ ] and unitary 
models chU[ ] is explained in the text below. These models 
exhibit a critical exponent close to 2/3. 





L/a 


In So 






-Fmin 


chO: RRH W/to = 5, zi 


= i.OO 






5 


12-96 


-1.60 ±0.41 


n Q1 Q _L A ni Q 

U.oio =t U.Uio 


i.iUU ± U.Uo 


n on 


5* 


12-96 


-2.3 ± 1.5 


U.oUi ± U.Uo 


1 1 v 4- n no 
i.i / ± U.Uz 


nil 
U.ii 


chO: RRH W/to = 6, zi 


1 TA 
= i. (\) 






5 


12-160 


-5.50 ± 0.20 


0.565 ± 0.006 


1.220 ± 0.003 


2.24 


5* 


24-160 


-4.17 ±2.02 


U. ( z4 ± U.U4D 


1 1 fin 1 nil 

i.ioU ± U.ii 


U.zo 


chO: RRH W/to = 7, zi 


= 2.23 






5 


12-128 


-3.65 ±0.18 


0.675 ± 0.005 


1.272 ± 0.002 


A SO 


5* 


12-128 


-4.4 ± 1.5 


0.693 ± 0.03 


1.204 ±0.02 


n 'HQ 


chO: RRH W/to = 8, zi 


= 2.65 






5* 


12-96 


-6.4 ±0.5 


0.628 ±0.01 


1.25 ±0.03 


1.71 


chS: RRH W/to = 4, 21 


= 1.97 






4 


8-160 


-5.4 ±0.1 


0.687 ± 0.049 


1.403 ±0.021 


7.40 


4* 


8-160 


-7.5 ±0.5 


0.692 ±0.015 


1.257 ± 0.004 


1.22 


4 


16-160 


-6.1 ±0.1 


0.676 ± 0.058 


1.409 ± 0.041 


6.40 


5 


16-160 


-3.4 ±0.1 


0.796 ± 0.078 


1.045 ± 0.033 


2.83 


5 


48-160 


-6.1 ±0.7 


0.664 ± 0.028 


1.186 ±0.025 


0.98 


chU[l/16, 5.0]: B = l/16/i/(2ea2), RRH W/to = 5, 












Z\ = 


= 2.47 


6 


32-96 


-1.74 ±0.37 


0.636 ±0.01 


1.27 ±0.13 


1.4 


7 


32-96 


-1.23 ±0.64 


0.648 ± 0.02 


1.17 ±0.17 


0.8 


chU[l/12, 0.5]: B = 1/12 h/(2ea'^), RMF / = 0.5 h/e 












21 = 


= 1.03 


4 


24-72 


-6.74 ±0.24 


0.567 ± 0.008 


1.75 ± 0.003 


3.504 


5 


24-72 


-4.14 ±0.65 


0.643 ± 0.02 


1.313 ±0.02 


2.142 


6 


24-72 


-3.45 ± 0.65 


0.663 ± 0.02 


1.332 ±0.23 


2.055 


7 


24-72 


-3.02 ±0.96 


0.676 ± 0.03 


1.411 ±0.33 


2.016 



2/3. Please note that the energies of the quantum Hall 
critical states belonging to the zcroth Landau band are 
splili^i by the disorder, and so they do not interfere with 
the chiral critical point at E' = studied here. 

To compare the new results with our previous data 
published in Ref. [l^ for the chiral unitary system with- 
out a constant magnetic field, we use the present fitting 
method also for the analysis of our previous data. In 
Table HI] we sec that the critical exponent for the model 
chU[0, 0.5] is consistent with the value 1/2 as reported 
previously;^ but differs significantly from the value 2/3 
obtained here for all other ensembles. Also, in order 
to check a possible disorder dependence of k, again a 
bricklayer model with zero constant magnetic field but 
with a stronger RMF disorder strength / = 0.7 h/e (la- 



TABLE II: The critical parameters for models with RMF dis- 
order and zero or small constant magnetic fields as obtained 
by the fitting procedure described in Section IIIII A critical 
exponent close to 1/2 is found. The energy interval used was 
\E\ < W~^°to. The labeling of the various models chU[ ] is 
explained in the text. 



Np 


L/a 


\nEo 






-Fmin 


chU[0, 0.5]: B = 0, RMF / = 0.5 h/e, zi = zl = 1.555 


6 


8-160 


-2.6 ± 1.6 


0.414 ±0.04 


2.49 ± 0.44 


0.91 


7 


8-160 


0.25 ± 2.5 


0.476 ± 0.05 


1.98 ±0.5 


1.43 


8 


8-160 


1.86 ± 1.1 


0.511 ±0.025 


2.01 ±0.44 


1.28 


chU[0, 0.7]: B = 0, RMF / = 0.7 h/e, zi = 1.22 


5 


24-96 


-4.5 ± 1.5 


0.515 ±0.04 


1.86 ±0.20 


2.43 


5 


32-96 


-3.5 ±1.5 


0.55 ±0.035 


1.68 ±0.11 


2.01 


5 


40-96 


-4.0 ± 1.7 


0.54 ±0.04 


1.65 ±0.10 


1.86 


chU[0, as 


: B = 0, RMF f/{h/e) = - 


-0.25,0.125, 


21 = 1.328 


5 


16-96 


-4.58 ± 1.8 


0.474 ± 0.04 


2.156 ±0.02 


1.36 


6 


16-96 


-1.66 ± 2.6 


0.535 ±0.05 


2.40 ±0.38 


1.22 


chU[l/10K, 0.5]: B = 


1/10000 h/{2ea^), RMF / = 


= 0.5 h/e. 










21 = 1.5559 


7 


24-96 


0.38 ± 4 


0.464 ± 0.07 


2.38 ± 1.14 


2.01 


8 


24-96 


1 ± 1 


0.475 ± 0.02 


3.16 ±0.75 


2.39 


chU[l/lK, 0.5]: B = 1/1000 /i/(2ea2 


), RMF / = 0.5 /i/e, 










21 = 1.557 


6 


24-96 


-0.7 ±3.7 


0.448 ± 0.063 


2.27 ±0.64 


1.82 


7 


24-96 


2.7 ±2.5 


0.511 ±0.043 


2.40 ±0.75 


2.52 


8 


24-96 


2.8 ± 1.4 


0.512 ±0.023 


2.73 ±0.44 


2.22 



beled chU[0, 0.7]) was calculated. We obtain the same 
value K ~ 1/2 and so confirming our previous results. In 
Fig. [T] we compare some of the raw numerical data for 
Zi{E,L) — Zi{E ~ 0,L), i = 1,2, with the fit given by 
Eq. ([9|). Some other models, listed in Tabic HIl which 
also exhibit an exponent k = 1/2 will be discussed in the 
next section. 



V. DISCUSSION AND CONCLUSIONS 

Our numerical data for two-dimensional models with 
various chiral symmetry confirm the logarithmic en- 
ergy dependence of the localization length reported 
previously 1^ Here, we find that this logarithmic depen- 
dence is determined by the critical exponent k, which is 
universal and close to 2/3 for chO, chS, and chU (with 
magnetic field) symmetry classes. The only exception is 
the chiral unitary system with zero constant magnetic 
field where k ~ 1/2. We have confirmed this previ- 
ous result by additional calculations for a RMF disorder 
strength / = 0.7 h/e. Therefore, the chiral unitary mod- 
els with and without constant magnetic field studied here 
exhibit different critical behaviors. 

One source for this outcome could be that the distri- 



5 



chU[0, 0.5] 
chU[0, 0.7] 




FIG. 1: (Color online) The difference -\z.,{E,L) - Zi{E = 
0, L) \ of the two smallest dimensionless scaling variables plot- 
ted as a function of the parameter x with corresponding fits 
(solid lines) for some systems listed in Table |T] and |TT] The 
previously obtained chU[0, 0.5] d&ta^ are shown only in part 
since x increases up to values of ~ 10 in this case. 
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FIG. 2: (Color online) The smallest scaling variable z\{E — 
0, /) versus RMF disorder strength /. The constant magnetic 
field is zero (red dots), B = 1/24 x h/{2ea'^) ( green squares), 
and B — 1/12 x h/{2ea?) (blue triangles), respectively. These 
results do not dependent on the system width in the range 
24 < L/a < 128 checked. 



bution of the random magnetic fluxes, vifhich was chosen 
to be symmetric about zero flux in the B = situation, 
turns into an asymmetric one due to the addition of a spa- 
tially constant magnetic flux in the finite magnetic field 
case. We have checked that an asymmetric distribution 
of RMF- fluxes still leads to a critical exponent 1/2 as 
long as the average magnetic flux remains zero. In these 
calculations, the random flux was taken to be —0.25 h/e 
with probability 1/3 and 0.125 ft-/e with probability 2/3 
so that the distribution of random fluxes is asymmetric 
albeit the mean value is zero. Using a 6-parametric fit as 
given by Eq. (10), we obtained that k — 0.535 ± 0.05 for 
16 < L/a < 96. The details for the antisymmetric RMF 
distribution chU[0, as] can be seen in Table llll 

Another reason for the different exponents could be 
that, for the largest system size studied, the disorder 
strength is still too small in the zero constant mag- 
netic field situation. Due to the periodicity, the effect 
of the random flux disorder — //2 < (px.y < //2 ap- 
pearing in the phases exp{~i9x^y+a-x,y) with Ox,y+a-x,y = 
'Px.y/h is bounded, being maximal 
for disorder strength / = 1.0 /i/e for the bricklayer lat- 
tice. The dependence of the smallest scaling variable 
zi{E = 0, /) on the disorder strength / > 0.01 is shown 
in Fig. @ for S = 0, p/q = 1/24, and p/q = 1/12, 
respectively. In the small disorder regime, the behav- 
ior is complicated and not important for the problem of 
different critical exponents discussed here. For large /, 
the zi{E = 0,/) decrease in all RMF models, indepen- 
dent of the presence or absence of a constant magnetic 
field, and become equal ~ 0.76 for / = h/e. Usually, a 



decreasing zi is equivalent to an increasing localization 
length A = 2Lj z\. A similar increase of the two-terminal 
conductance with increasing RMF-disorder strength has 
been reported previously)^ The behavior of the RMF- 
models is opposite to the chO and chS situation where 
z\(E = 0, W^) becomes larger with increasing disorder 
strength W . However, when comparing the various val- 
ues of zi{E = 0) obtained for different models (see Ta- 
ble |l|), which are related to the effective disorder strength, 
the above suggestion that too weak a disorder was applied 
in the RMF-model without B-field seems to be unlikely, 
and we do not believe that this is the correct explanation 
for the observed difference in the critical exponents. 

Another distinction between the two chU cases is the 
behavior of the density of states p{E) very close to 
E = Q. Although both the chU models with and with- 
out strong magnetic field exhibit a narrow depression^ 
with p{E) going to zero at the critical point, the disorder 
dependence is quite different. With increasing random- 
magnetic-flux disorder strength, the energy range of the 
depression gets narrower in the absence of a constant 
B while it increases when a strong perpendicular mag- 
netic field is present. The appearance of such 'micro- 
gaps' has been attributed to the non-perturbative ergodic 
regimeii^i^ of the chiral model, i.e., the localization 
length greatly exceeds the sample size so that for time 
scales large compared to the diffusion time the particle 
extends over the whole sample. A similar behavior was 
observed also for the chS model.— We note, however, that 
DOS depressions ('microgaps') are absent in our numer- 
ical calculations performed for the chO case. For strong 
disorder, a narrow peak appears instead ai E = 0^ pre- 
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sumably originating from the different type of disorder 
(real random flopping) that in contrast to the complex 
random phases in the RMF models can lead to isolated 
sites having eigenenergies around E = 0. This extra peak 
is observed for both the bricklayer and the square lattice 
systems. Nevertheless, it is unclear whether the different 
disorder dependence of the DOS around the critical point 
can account for the observed difference in the critical ex- 
ponents. 

Finally, we would like to remind the reader that the 
presence of a strong magnetic field causes a special topo- 
logical term in the appropriate field theories i^i^i^ This 
extra term is responsible for the occurrence of current 
carrying states that are essential for the explanation of 
the integer quantum Hall effect, and so it accounts for the 
difference between the ordinary Gaussian unitary ensem- 
ble and the quantum Hall situation, i.e., the appearance 
of critical electronic states in the latter case. 

Furthermore, the strong magnetic field may effectively 
restore the particle-hole symmetry that is usually lack- 
ing when time-reversal symmetry is broken. The recov- 
ery of the particle-hole symmetry can take place when 
the scattering between electronic states belonging to dif- 
ferent disorder broadened Landau bands is suppressed 
by a large energy separation. Thus, a single Landau 
band model containing a chiral quantum critical point 
at i? = is effectively created in the chiral unitary case 
with an additional very strong, spatially constant, per- 
pendicular magnetic field. We therefore suggest that the 
presence or absence of particle-hole symmetry is respon- 
sible for the different critical exponents found in our nu- 
merical calculations of various chiral unitary models. To 
check this view, we investigated also chiral unitary mod- 
els with weaker constant magnetic fields p/q = f/fOOOO 



and p/q = 1/1000 (see chU[l/10K, 0.5] and chU[l/lK, 
0.5] in Table HI)) and found critical exponents ~ 1/2 as in 
the magnetic field free case. 

In conclusion, we calculated the two smallest Lya- 
punov exponents of lattice models belonging to the chi- 
ral orthogonal, chiral symplectic, and to the chiral uni- 
tary class in the presence of a strong spatially constant 
magnetic field. We found in all cases a critical expo- 
nent K ~ 2/3 that governs the divergence of the local- 
ization length at the band center E = 0. This result 
is in agreement with previous calculations for the di- 
verging energy dependence of the density of states.— 
It is, however, at variance with the earlier analytically 
obtained results.— i^i While reasons for this difference 
were suggested in Refs. ( ll l20ll2^l23 ) . the origin of the 
numerically^^ obtained value k ~ 1/2 for chiral unitary 
lattice models in the absence of a strong perpendicular 
magnetic field is still an open question. The latter result 
was again corroborated in the present study and sug- 
gested to be due to the absence of particle-hole symmetry 
in the RMF-model without a strong spatially constant 
magnetic field. However, due to the limited achievable 
system size in the numerical calculations, a magnetic 
field induced transition from the ballistic to the diffu- 
sive regime as discussed above, cannot be completely ex- 
cluded. 
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